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Abstract  of  Dissertation  Presented  to  the  Graduate  Council 
in  Partial  Fulfillment  of  the  Requirements  for  the 
Degree  of  Doctor  of  Philosophy 

A SPECTROMETRIC  INVESTIGATION  OF  AN 
ELECTRODELESS  DISCHARGE 

By 

Kynric  Martin  Pell 
August,  1967 

Chairman:  Dr.  M.  H.  Clarkson 

Major  Department:  Aerospace  Engineering 

This  dissertation  is  concerned  with  the  experimental 
determination  of  the  electron  temperature  distribution  in 
a nonequilibrium,  low  pressure,  electrodeless  discharge. 

In  view  of  the  lack  of  excitation  function  data  for  com- 
plex atoms,  such  as  argon  which  was  used  in  this  study,  an 
approximate  semi-corona  model  was  used.  The  temperature 
dependence  of  excited  state  population  densities  for 
neutral  atoms  and  singly  ionized  ions  predicted  by  the 
approximate  model  was  used  to  obtain  a relationship  between 
the  intensity  ratio  of  two  emission  lines  and  the  temper- 
ature . 

Integrated  intensity  distribution  data  were  inverted 
numerically  to  obtain  the  radial  dependence  of  the  in- 
tensity . 

Ratios  of  the  intensity  of  an  ion  emission  line 
and  an  atom  emission  line  at  25  radial  locations  in  the 
plasma  were  used  to  obtain  the  temperature  distributions. 


xi 


Maximum  temperatures,  found  on  the  axis  of  symmetry, 
ranged  from  1.24  to  1.54  ev.  for  a range  of  operating 
pressures  from  50  to  1000  microns  of  mercury.  The  radial 
dependence  of  the  temperature  was  found  to  be  quite  weak. 
The  temperature  decreases  approximately  15%  over  a region 
equal  to  2/3  of  the  radius  of  the  plasma. 
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CHAPTER  I 


INTRODUCTION 

An  electrodeless  discharge  is  a plasma  which  is 
formed  by  the  excitation  of  a gas  by  either  c.w.  or  pulsed 
radio  frequency  (r.f.)  fields.  The  discovery  of  the  elec- 
trodeless discharge  is  usually  attributed  to  Hittorf,^  in 
the  year  1884.  For  over  40  years  following  its  dis- 
covery, a controversy  raged  regarding  the  nature  of  the 

1 2 

discharge.  Hittorf  and  J.  J.  Thomson  contended  that  the 

discharge  was  due  to  electromagnetic  induction.  Townsend 
3 4 

and  Donaldson,  as  well  as  Lehrmann,  among  others,  main- 
tained that  the  discharge  was  due  to  the  large  potential 

difference  between  the  ends  of  the  current  carrying  coil. 

5 

In  1929,  K.  A.  MacKinnon  described  two  different  modes  of 
operation  he  found  in  the  electrodeless  discharge  which 
were  primarily  dependent  upon  working  gas  pressure,  coil 
configuration,  and  coil  current.  He  found  that  Thomson's 
work  was  confined  to  a small  range  of  pressures  involving 
ring- type  discharges  exclusively.  MacKinnon  concluded 
that  ring  discharges  were,  in  fact,  due  to  electromagnetic 
fields.  On  the  other  hand  Townsend  and  Donaldson  treated 
a glow  discharge  over  a wide  range  of  pressures  and  ap- 
parently failed  to  realize  the  difference  between  the  ring 
and  glow  modes  of  operation.  MacKinnon  concluded  that  the 
glow  mode  of  operation  was  either  largely,  or  entirely, 
due  to  electrostatic  excitation. 
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This  resolved  the  controversy  and  it  was  1960  before 

0 

Reed  applied  electromagnetic  induction  to  a flowing  gas 

at  atmospheric  pressure  to  create  an  electrodeless  plasma 

torch.  Quantitative  temperature  determinations  made  by 

Reed,  using  spectrographic  techniques  and  the  Fowler-Milne 

method,  resulted  in  the  temperature  distribution  shown  in 

7 

Figure  1.  In  1963  Reed  investigated  static  plasmas  of 
several  gases  over  a range  of  pressures  from  .005  mm  of 
mercury  to  atmospheric  pressure  using  spectroscopic  tech- 
niques. He  observed  two  different  modes  of  operation  which 
were  reflected  in  different  spectra  obtained  at  various 
pressures  and  plate  powers.  However  he  did  not  make  fur- 

g 

ther  quantitative  temperature  predictions.  Pell,  in 

1963,  utilized  a stagnation  point  heat  transfer  probe  to 

measure  the  temperature  distribution  in  an  atmospheric 

pressure  plasma  torch  yielding  a significantly  different 

temperature  distribution  than  obtained  by  Reed.  Spectro- 

9 10 

scopic  work  by  Gold'farb  and  Dresvin,  Pettit,  and 
Johnston^  among  others  has  confirmed  the  lower  temper- 
atures obtained  by  Pell  for  the  atmospheric  pressure 
plasma  torch.  A typical  temperature  distribution  for  such 
a torch  is  shown  in  Figure  1. 

At  pressures  below  atmospheric,  probe  studies  have 

12  13 

been  made  by  Clarkson,  et  al.,  ' and  Keefer,  et  al. 

14 

Some  spectrometric  measurements  were  made  by  Porter 

wherein  a corona  model  of  the  plasma  was  used  in  the  in- 

15  16 

terpretation  of  the  data.  Sovie  ’ has  applied  both 
probe  and  spectrometric  techniques  to  a flowing  helium 
plasma  in  the  pressure  range  of  0 to  100  microns  obtaining 
reasonable  agreement  between  the  two  techniques. 

The  problem  associated  with  any  spectrometric  diagnostic 
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attempt  on  the  low  pressure  discharge  is  that  it  is  diffi- 
cult to  determine  the  number  density  of  atoms  or  ions  in 
the  various  excited  states.  Sovie's  success  with  helium  was 
dependent  on  the  hydrogenic  approximation  and  the  fact  that 
excitation  cross-sections  for  such  ions  have  been  calcu- 
lated and  experimentally  verified. 

The  problem  addressed  in  this  work  is  the  temperature 
determination  in  a low  pressure  (50  to  1000  microns)  argon 
electrodeless  discharge.  A semi-corona  model  of  the  plas- 
ma is  used  in  this  study.  A functionally  simple  semi- 

17 

corona  ionization  equation  described  by  Oertel  is  applied 
to  an  argon  electrodeless  discharge.  Non-equilibrium, 
steady-state  populations  of  the  various  excited  states  of 
the  atom  and  ion  predicted  by  Oertel's  method  are  used  to 
obtain  a relationship  between  the  intensity  ratio  of  two 
emission  lines  and  the  electron  temperature. 

Integrated  spectrometric  intensity  data  for  four  argon 
emission  lines  are  inverted  using  Abel's  equation  to  ob- 
tain radial  intensity  distributions  in  the  discharge. 

These  radial  distributions  are  used  to  obtain  local  values 

°f  I7503/I4806  Md  I7514/I4347  fr°m  WhlCh  temperature 
distributions  are  obtained. 

Temperatures  were  obtained  for  a range  of  input  power 
levels  from  220  to  740  watts  in  the  pressure  range  from  50 
to  1000  microns.  These  results  are  presented  in  Table  I. 


CHAPTER  II 


PLASMA  MODELS 


Discussion 

The  basic  problem  associated  with  any  spectroscopic 
technique  is  relating  the  intensity  of  radiation  to  the 
plasma  properties. 

There  are  currently  three  microscopic  models  of  a 
plasma  in  use  for  the  interpretation  of  spectrometric 
data.  These  models  which  are  termed  thermal,  corona, 
and  semi-corona  (or  collisional-radiative)  are  assumed 
to  apply  to  plasmas  in  different  states.  The  rate  pro- 
cesses for  populating  and  depopulating  the  various  excited 
states  in  the  plasma  are  proportional  to  the  number  density 

of  electrons  N and  the  electron  temperature  T . Popu- 
® e 

lation  by  photon  induced  excitations  or  ionizations  is 
usually  neglected  on  the  assumption,  or  experimental 
verification,  that  the  plasma  is  optically  thin.  The 
domains  of  applicability  of  the  various  models  can  be 
represented  on  a graph  of  N vs  T for  any  particular 

6 G 

type  of  gas  and  this  has  been  done  in  Figure  2 for  argon. 

The  criteria  used  in  constructing  such  a graph  are  pre- 
sented in  the  following  sub-sections.  A description  of 
each  of  the  three  models  follows. 

Thermal  Plasma 

Thermal  plasmas  are  those  in  which  the  population  of 
excited  states  is  in  accord  with  the  Boltzmann  distribution. 
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That  is 

Hz  = 3f>  exp[-Ep/kT\ 

N uTT) 

(2.1) 

where 

U(T)=  Y.3P  e*P  [-  er/k  r) 

p (2.2) 

N = total  number  density  of  atoms  in 

all  excited  states  and  the  ground  state. 

Np  = number  density  of  atoms  in  excited  state  p. 
g = statistical  weight  of  state  p. 

r' 

Ep  = excitation  energy  of  state  p. 

k = Boltzmann  constant. 

T = Temperature. 

U(T)  = Partition  function. 

For  complete  thermodynamic  equilibrium  to  exist,  detailed 
balancing  is  required.  This  condition  is  rarely  found  in 
laboratory  plasmas.  However,  the  concept  of  local  thermo- 
dynamic equilibrium  (LTE)  may  be  applicable  when  detailed 
balancing  is  not  satisfied. 

Whenever  local  thermodynamic  equilibrium  exists,  the 
densities  in  specific  quantum  states  are  those  existing 
in  a system  in  complete  thermodynamic  equilibrium  which 
has  the  same  mass  (total)  density,  temperature,  and 
chemical  composition  as  the  LTE  system.  The  temperature 
referred  to  is  that  of  the  specie  dominating  the  reaction 
rates.  In  general  this  is  the  electron. 
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The  distribution  of  excited  states  in  LTE  is  given 
by 

Nf>  m 9p  [~Tp/kTe\ 

N U(Te) 

(2.3) 

where  is  the  electron  temperature.  It  is  not  necessary 

for  the  ions,  nor  the  atoms,  to  have  the  same  temperature 

as  the  electrons.  However,  T > T and  T . 

e " g i 

In  an  LTE  plasma  ionization  equilibrium  is  expressed 
by  the  Saha  equation 


NZ**N  e z UZfl 

kTe  ‘ 

34 

exp 

<vz  *.*  u1 

kTe 

where 


(2.4) 


N 

N 

e 

00o 

eh 

AEE 


number  density  of  ions  in  ionization 
stage  Z.  (Z  = 1 for  atoms) 

number  density  of  electrons. 

first  Bohr  radius, 
ionization  energy  of  hydrogen. 


reduced  ionization  potential. 

The  hydrostatic  pressure  of  the  gas  is  given  by  Dalton's 
law  of  partial  pressures 


p = 

1 (2.5) 

where  the  i includes  all  species  and  stages  of  ionization. 
The  additional  requirement  of  quasi-neutrality 
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He  = X.Ui-liNt 

(2.6) 

where  i refers  to  stage  of  ionization,  completes  a set  of 
equations  which  can  be  used  to  determine  the  composition  of 
an  LTE  gas . 

For  this  model  to  apply  to  a gas,  the  electron  density 
must  be  sufficiently  large  to  insure  that  the  population 
of  all  states  is  dictated  by  electron  collisional  excita- 
tion and  de-excitation  rates.  Radiative  losses  are  assumed 
to  be  negligibly  small. 

Corona  Plasma 

The  most  extreme  departure  from  LTE  is  described  by 
the  corona  model  which  has  been  proposed  as  a model  of  the 
solar  corona.  (See,  e.g.  reference  18).  In  LTE  the 
process  of  collisional  ionization  is  balanced  by  three- 
body  recombination  , 


e t Nz 


. z*‘1 

N + e t e 


(2.7) 

As  the  number  density  of  electrons  decreases,  the  process 
of  radiative  recombination  becomes  large  compared  to  the 
three-body  recombination  process,  i.e. 


e t"  A/*'1  — N*  + hv 

(2.8) 

In  this  model,  it  is  assumed  that  the  electrons  have 
a Maxwellian  distribution  of  velocities  and  that  the  mean 
energy  of  the  ions  and  atoms  is  equal  to,  or  less  than, 
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that  of  the  electrons.  In  addition  it  is  assumed  that  the 
population  of  all  excited  states  is  negligible  compared  to 
the  ground  state.  An  equation  relating  the  collisional 
ionization  rate  to  the  radiative  recombination  rate  may 
now  be  written  as 

Ne  Mj  S(Te,2,3)-He  Nafl<z(Te,2ti,q) 

J (2.9) 

or 


Na  m cl  ( Te  jZ+1  ,g) 

ure)z)9) 

(2.10) 

Results  of  detailed  calculations  of  S by  Burgess, 

which  were  based  on  the  Coulomb-Born  approximation,  are 

18 

given  by  McWhirter.  Also,  a general  formula  for  0C  is 

18 

given  by  McWhirter.  The  expressions  for  S and  (£  may  be 
used  to  write  an  ionization  balance  equation  of  the  form 


N?1  $Te3/+  VhTe. 

<2.11 ) 

where  £ is  the  number  of  equivalent  electrons. 

The  population  of  the  pth  excited  state  may  be  re- 
lated to  the  ground  state  by  the  expression 

A*/  = A'c  XJZ_, 

n3  f\  (z,p) 


(2.12) 
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where  X(Z,g,p)  is  the  collision  excitation  function  for 
the  temperature  and  A(Z,p),  sometimes  written  as  Ap, 
is  the  total  transition  probability  for  the  state  p. 

Application  of  this  model  requires  that  the  electron 
number  density  be  sufficiently  low  so  that  the  collisional 
de-excitation  term  can  be  neglected  in  the  ionization 
balance  equation  (2.9). 

Collision-Radiative  Plasma  (Semi-Corona) 

When  the  electron  density  is  too  low  to  insure  LTE 
and  too  high  to  apply  the  corona  model,  the  electronic 
state  of  the  atoms  and  ions  in  the  plasma  must  be  deter- 
mined from  a system  of  rate  equations.  That  is,  the  steady 
state  population  of  an  excited  level  is  a function  of  both 
radiative  and  collisional  processes.  The  assumptions 

required  to  set  up  such  a system  of  equations  are  listed 

. , 18 
below : 

(1)  The  free  electrons  must  have  a Maxwellian 
distribution 

(2)  Ionization  is  by  electron  collision  from  the 
bound  levels  of  the  lower  stage  of  ionization 

, and  is  only  partially  balanced  by  three-body 

recombination 


A/Z  . A ,2il 

N P te  ^ N j + e + e 

(3)  Transitions  between  bound  levels  may  be 
induced  by  electron  collisions 


(2.13) 


f/jo  f e ** — N g + & 

(4)  Radiation  is  emitted  when  an  electron  in  a 
higher  bound  state  p,  spontaneously  decays 
to  a lower  state  g and  when  a free  electron 
makes  a transition  to  a bound  level 
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(2.15) 


(2.16) 


Using  the  above  assumptions  the  equations  describing 
the  rates  at  which  the  various  excited  states  are  populated 
may  be  written 


where  ^(T^,  Z+l,  p)  is  the  collisional  recombination 
coefficient . 

In  general  the  solution  of  a set  of  equations,  typi- 
fied by  (2.17)  for  all  states  of  the  atoms  and  ions  in  a 
plasma  is  not  possible  even  when  the  number  of  states  is 
reduced  to  a finite  number  by  the  reduced  ionization  po- 
tential. The  problem  is  not  the  number  of  equations  but 
the  lack  of  rate  coefficients  S(Te,Z,p),  A(Z,p,g),  etc. 
for  complex  atoms  or  ions.  Detailed  calculations  for 


(2.17) 
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hydrogen  and  hydrogen-like  ions  have  been  made  by  several 
17 

researchers  using  several  assumptions  in  addition  to 

the  four  given  above. 

17 

Oertel  has  proposed  an  approximation  for  the  semi- 
corona model  which  assumes  that  the  electronic  excitation 
of  the  atoms  or  ions  can  be  described  by  two  different 
Boltzmann  distributions  while  the  ionization  equilibrium 
is  given  by  a corrected  Saha  equation. 

If  the  electron  number  density  in  an  arbitrary  plasma 
is  reduced,  departures  from  LTE  will  start  at  the  lowest 
energy  levels  of  the  atom  or  ion.  For  any  steady  state 
in  the  collision-radiative  or  semi-corona  domain,  the 
ratio  of  population  densities  of  the  higher  states  are 
given  by  the  Boltzmann  relation 


Nm 


exp 


-(£%  - Em) 
kTe  . 


(2.18) 

That  is,  the  higher  states  are  in  LTE.  For  the  lower  states, 

a correction  factor  T is  introduced  so  that  (2.18)  be- 

nm 

comes 


A/f 

Nm 


(2.19) 

2"nm  is  obtained  by  requiring  (2.19)  to  predict  the  same 

ratio  as  the  corona  model.  An  expression  for  "T  based  on 

u nm 

17 

hydrogen  and  hydrogenic  ions  is  obtained  by  Oertel  . 

However,  a simpler  form  is  given  by  Schneider  and  Oertel^ 
which  shows  better  agreement  with  detailed  calculations 
by  Bates,  et  al.20 
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This  expression  for  T is 

nm 


r,3  Ve  Iff 

Eh  ‘ 

3 r^j 

1 

8 Nr 
J 

L£l. 

i Vi 


J ftr  T'<.1 


'Fntn  - i.  } -for  ^ j_ 


(2.20) 

where  <£  = fine  structure  constant,  and  £0  = collision  limit 
to  be  defined  later . The  modified  Saha  equation  used  in 
this  approximation  is 


Nz  CC*  U2 

o 


'-(eI  ZA£j[ 
kTe  . 


(2.21) 

The  collision  limit  is  a state  with  principal  quantum 
number  n’  where  the  excitation  level  of  n'  is 


E(  n')  = £ „ - £0 

(2.22) 

In  general  n'  is  nonintegral.  However,  there  is  no 
necessity  to  associate  n'  with  a particular  energy  level 
of  the  atom  or  ion  since  one  only  desires  to  know  the 
regions  of  applicability  of  the  various  Boltzmann  equations 
The  state  n*  is  determined  by  equating  the  collisional 
excitation  rates  to  the  radiative  decay  rates  and  solving 
for  the  energy  at  which  this  situation  will  occur.  The 
form  of  <£0,  as  defined  in  (2.22),  is  given  by  Schneider 


19 

and  Oertel  as 

= Eh 


3 A/e  (CC„J3  f Eh  jz  Z 


"l  2/a 


Eli 


(2.23) 


13 


„ • 21 
Gnem  as 


The  reduced  ionization  potential  A E^,  is  given  by 


aEZ  = z(z+l)Ei 


■H 


Sir  Eh  (Ut+  Zz-  tJ;) 

jfr^  j J. 


(2.24) 


However,  a much  simpler  form  proposed  by  Schneider22  has 
been  used  here,  namely 


aE 


« - o.&Zx/o  ^(N  ) +0'37X/0? ’[Me}/z 

e 1 TJ  c 


2.25) 

Regions  of  Validity  for  Local  Thermodynamic  Equilibrium 

In  order  that  a stationary  homogeneous  plasma  be  in 
local  thermodynamic  equilibrium,  a minimum  electron  density 
requirement  must  be  imposed  to  insure  that  radiative  tran- 
sitions between  two  states  are  less  than  10%  as  frequent  as 
de-excitations  by  electron  collisions.  Since  radiative 
rates  tend  to  decrease  with  higher  levels  of  excitation 
and  whereas  collisional  rates  tend  to  increase,  LTE  is 
more  probable  for  higher  excited  states  if  the  electron 

number  density  is  too  low  to  permit  complete  LTE.  Schneider 
19 

and  Oertel  give  a minimum  electron  density  to  insure 
LTE  between  excited  states  m,  located  above  a limiting 
state  n',  and  the  continuum  for  a non-hydrogenic  atom  in 
the  form 


He.  > / ?r/0 
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This  condition  was  used  to  set  the  lower  boundary  of  the 
LTE  region  on  Figure  2.  If  the  plasma  is  inhomogeneous, 
then  the  mean  free  path  of  atoms  between  equilibrating 
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collisions  must  be  short  compared  to  the  distance  over 
which  the  plasma  conditions  change  appreciably.  Further 
discussion  of  this  criterion  may  be  found  in  reference  21. 

Regions  of  Validity  for  Corona  Equilibrium 

The  corona  domain  may  be  defined  as  that  domain 
where  the  rate  coefficients  OC t and  St  for  the  thermal 
(higher  excitation)  levels  are  much  less  than  the  cor- 
responding coefficients  (£,  and  S for  the  continuum.  This 

23 

definition  has  been  used  by  Wilson  ' to  establish  an 
upper  limit  on  the  electron  density  given  by 

(ATJ* 

This  is  the  condition  used  in  the  graph  of  Figure  2 in 
order  to  establish  the  corona  limit. 

Regions  of  Validity  for  Collisiona 1-Radiative  or 
Semi-Corona  Domain 

As  may  be  seen  from  Figure  2,  the  region  between  the 
two  limits  defined  above  is  the  collisional-radiative 
domain . 


CHAPTER  III 


THE  MODEL  OF  THE  ARGON  PLASMA 


Probe  data  taken  by  Keefer,  et  al.,  indicate  that 
the  electrodeless  discharge  in  argon  operates  at  an 
electron  temperature  in  excess  of  1.0  volt  with  a maxi- 
mum electron  number  density  on  the  axis  of  approximately 
12  -3 

5.0  X 10  cm  . In  view  of  the  discussion  in  Chapter 
II,  it  appears  that  the  appropriate  microscopic  model  for 
the  low  pressure  electrodeless  discharge  is  the  semi- 
corona model.  The  lack  of  either  experimental  or  theo- 
retical cross-sections,  excitation  functions,  and  other 
atomic  parameters  for  complex  atoms,  such  as  argon,  make 
the  solution  of  the  rate  equations  a practical  impos- 
sibility. Two  approximate  forms  for  the  semi-corona 
region  are  available.  The  approach  due  to  Oertel  de- 
scribed in  Chapter  II  is  used  in  the  following  analysis 

because  of  its  simplicity.  A considerably  more  detailed 

21 

approach  is  given  by  Griem.  The  development  of  Oertel’s 
approach  for  an  argon  plasma  follows. 

The  intensity  ratio  of  two  emission  lines  in  a 
steady  state  plasma  is  given  by  the  equation 


(3.1) 
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. _(1)  . 

where  I is  the  intensity  of  spectral  line  (1)  having 
frequency  and  originating  from  transitions  between 
level  q and  level  p.  A is  the  transition  probability 
associated  with  the  levels  p and  q.  The  integration  is 
over  the  extent  of  the  plasma  along  the  line  of  sight  of 
the  detector.  If  the  experimental  data  can  be  unfolded 
so  that  local  intensity  ratios  can  be  obtained,  the  in- 
tegration may  be  eliminated  and  equation  (3.1)  may  be 
written 


Iw  v,  A/"1 
r«>  = W*Z  AW 

(3.2) 

where  the  intensity  ratio  is  a function  of  radial  location. 

The  problem  of  unfolding  the  intensity  data  is  discussed 

in  Chapter  IV.  In  general,  the  ratio  N^V  is  a 

Q n 

strong  function  of  the  temperature  if  the  excited  states 
q and  n are  in  different  ionization  stages.  Using  a semi- 
corona model,  the  ratio  will  be  related  to  the 

temperature  so  that  an  intensity  ratio  measurement  can 
be  used  to  determine  the  temperature. 

The  number  density  in  a particular  excited  state 
n,  of  ionization  stage  (Z),  can  be  related  to  the  total 
number  density  of  the  element  by  the  equation 


A/'n1  = //J,*’  tfU> A/ 

N,z>  V 


(3.3) 
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For  a system  in  equilibrium 


Nn 


N 


(z) 


Ni 


(2) 


(Z) 


N™  U,v(kTe) 


(3.4) 

where  g refers  to  the  ground  state.  The  usual  assumption 
made  in  analyses  of  this  type  is  that  the  partition  func- 
tion may  be  replaced  by  its  first  term,  so  that 


Mi 
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Nn  _ 

/v,z>  Uf 


(3.5) 

To  take  into  account  nonequilibrium  populations  a cor- 
rection factor  ng  is  introduced  into  the  Boltzmann 
relation  which  gives 
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(3.6) 


Equation  (3.3)  may  now  be  written  in  the  form 


(Z) 

Nn  = 


T»3  eXP 
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(3.7) 


(Z+l)  (Z) 

The  ratio  Nv  / N may  be  written  in  the  form 
n P 


(Z+l)  (2tl)  (ZJ 

Nn  yn  Jf  exp 

~r~'  - - * (ZWJz)  ‘ 
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( Z+l ) (Z ) 

The  ratio  of  N / Nv  is  given  by  the  corrected 
Saha  equation  (2.21).  Figure  3 shows  the  dependence  of 
"V  ^ on  temperature  for  argon  at  electron  densities 


.10 


.12 


.14 


-3 


of  10  , 10*“,  and  10~"  cm  . The  values  of  the  Saha 

correction  factor,  T , and  the  collision  limit,  £0  , re- 
quired to  construct  Figure  3 were  determined  using  equa- 
tions (2.20)  and  (2.23).  The  results  are  presented  in 
Figures  4 and  5,  respectively.  All  of  the  emission  lines 
used  in  this  study  originated  from  levels  lying  below  the 
collision  limit,  n'.  Therefore,  the  form  of  the  Boltz- 
mann equation  used  in  determining  the  populations  of  the 
appropriate  upper  states  of  the  atoms  and  ions  includes 
the  correction  factor  ^ . For  the  atoms,  this  correction 
factor  is  the  same  as  that  previously  used  in  the  Saha 
equation.  A similar  graph  was  made  for  the  ions.  How- 
ever, it  is  not  presented  here.  Using  equation  (3.8), 


the  correction  factor 


r <2> 


ng 

, and  the  ratio  N^2V  N^,  a graph  of  N N ^2V 

PQ.  n p 

where  n corresponds  to  the  107054.32  cm-  level  of  neutral 
argon,  (Arl),  and  p corresponds  to  the  157234.93  cm-1 
level  of  singly  ionized  argon,  (Aril),  is  plotted  for 
various  electron  densities  on  Figure  7.  A similar  cal- 
culation for  the  Arl  level  108722.67  cm-1  and  the  Aril 
level  155044.07  cm  1 gives  the  same  result.  Thus,  Fig- 
ure 7 represents  the  ratio  N N ^2^for  the  two  popu- 

n p 

lation  ratios  used  in  this  study. 

Using  equation  (3.2),  the  intensity  ratios  I__rt_,I.or._ 

7 duo/  4oUb 

and  I7514'/l4347  were  plotted  as  a function  of  temperature 

in  Figure  8 for  an  assumed  number  density  of  lO^2  elec- 

-3  4 -3 

trons  cm  . A variation  in  electron  density  of  10  cm 


the  correction  factor 
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results  in  a temperature  change  of  .15  ev.  or  approxi- 
mately 10%  for  the  temperatures  determined  in  this 

investigation.  The  transition  probabilities  used  in  this 

24 

study  were  those  given  by  Olsen. 


CHAPTER  IV 


RADIATION  FROM  A FINITE  SOURCE 
External  Optical  System 


To  make  relative  or  absolute  intensity  measure- 
ments on  an  inhomogeneous  finite  radiating  source,  an 
external  optical  system  is  usually  used  to  image  the 
source  on  the  entrance  slit  of  a spectrometer.  Since  the 
source  is  three  dimensional,  the  question  arises  as  to 
where  the  optical  system  should  focus  in  the  source  and 
indeed  how  does  one  interpret  the  data  since  they  are 
obtained  from  an  integrated  intensity  over  the  source 
in  the  field  of  view  of  the  optical  system.  An  analyti- 
cal treatment  of  the  radiation  from  a finite  source  is 
considered  in  the  following  section  which  is  then  fol- 
lowed by  a description  of  an  experimental  verification 
of  the  results. 

Irradiance  of  the  Image  of  an  Extended  Source 

Consider  the  system  shown  in  Figure  9 where  the 
source  is  of  radiance  N and  its  image  is  of  radiance  N' . 

We  wish  to  find  the  irradiance  H at  any  point  in  the  image 
space.  Let  d^  be  an  element  of  area  at  a distance  X from 
the  source.  Its  image  is  d d Consider  a pencil  of  rays 
of  radiant  intensity,  J,  included  in  the  angle  d 6 at 
an  angle  & from  the  optical  axis.  If  Lambert's  cosine 
law  is  obeyed,  then  the  irradiance  dH  at  dd  is  given  by 


20 


21 


dH  = dj_  cos  e 
X 2 sec * 6 

But  dJ  = Zft r dr  H cos  e 


so  that 


dH  - zirrdr // cos2,e 
Xz  sec 2 e 


2 

Note  that  r = X tan  0 and  that  dr  = X sec 

r 6 

Therefore  H = / ZffH Sib/  & CoSB  d& 

J o 


& • 


d e 


= 7X  M s/ye0 

The  Abbe  sine  condition  gives  the  magnification  as 

/r?  = JL  - U S/A/6 

1 JU's/A/0' 

Squaring  and  rearranging  gives 


S//7Ze  - mz 3//vZ&' 

Conservation  of  energy  in  the  optical  system  yields 


H - rrjZH  ' 

Substituting  (4.4)  and  (4.5)  into  (4.3)  gives 


/rj  *//  = vA/mz{ M J ^s/a/^g' 


(4.1) 


(4.2) 


(4.3) 


(4.4) 


(4.5) 


(4.6) 
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For  source  and  receiver  planes  both  in  air,  M -M 
and  the  Law  of  Irradiance  may  be  written 

//  = 7r  /V  8/nZ  &'  (4  7) 

This  result  is  independent  of  the  location  of  d^  ’ 
so  that  if  we  assume  that  this  plane  coincides  with  the 
plane  of  the  entrance  slit  we  see  that  the  irradiance  of 
the  entrance  slit  by  any  surface  in  the  source  is  given 
by 

N ’=  7rN  s/n2 8 ' 

(4.8) 

If  it  is  not  desired  to  resolve  the  intensity  dis- 
tribution spatially,  it  is  not  necessary  to  use  an  imaging 
optical  system  for  the  external  optics. 

Experimental  Verification  of  the  Law  of 
Irradiance  for  an  Extended  Source 


A glow  discharge  tube  was  mounted  on  a movable  stand 
which  was  attached  to  an  optical  bench.  The  discharge 
tube  could  be  moved  up  and  down  (scanned  across  its  di- 
ameter) in  .01  mm  increments.  An  image  of  the  discharge 
tube  was  focused  on  a slit  by  a single  lens.  A photomul- 
tiplier located  behind  the  slit  recorded  the  total  inten- 
sity of  radiation  from  the  field  of  view  of  the  optical 
system.  An  adjustable  stop  was  used  to  vary  the  f number 
of  the  optical  system.  For  a range  of  / numbers  from  9 
to  36  and  an  out  of  focus  condition  of  - 10%  of  twice  the 
focal  length  of  the  lens  the  intensity  distribution  for  the 
discharge  tube  as  recorded  by  the  photomultiplier  remained 
essentially  unchanged.  The  deviations  in  the  intensity 
distributions  for  the  various  f numbers  and  focusing 
conditions  was  less  than  1%  as  shown  in  Figure  10. 
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Spatial  Resolution  of  Cy lindrica lly 
Symmetric  Sources 


It  may  be  seen  from  the  law  of  irradiance  that  the 
amount  of  light  entering  the  entrance  slit  of  the  spec- 
trometer from  an  extended  uniform  source  is  independent 
of  the  position  of  the  source  on  the  optical  axis  of 
the  system.  If  the  source  is  in  focus  within  a reason- 
able experimental  error  and  the  plasma  properties  can  be 
considered  constant  across  a plane  perpendicular  to  the 
optical  axis  which  is  bounded  by  the  field  of  view  of 
the  optical  system,  the  irradiance  of  the  entrance  slit 
is  given  by 


H = rr  s/nz & '/  //M  dx 

^ A 


'0  (4.9) 

where  the  integration  is  along  the  line  of  sight  (optical 

axis)  through  the  plasma.  For  spectroscopic  work  one  is 

interested  in  relating  spectral  irradiance,  that  is  power 

-2 

per  unit  wavelength  interval  per  cm  incident  on  the 
entrance  slit,  to  the  spectral  radiance  of  the  source, 


i . e . 


ft H - ft  $jnz&'  f Nh  (X)  JX 


'o  - (4.10) 

The  quantity  measured  by  the  spectrometer  is  H*  K , 

/l  i 

where  is  an  experimental  constant  which  depends  on  the 
system  response  of  the  spectrometer  and  the  geometry  of 
the  external  optical  system.  Assuming  the  geometry  of  the 
system  is  held  constant,  only  the  wavelength  dependence 
of  K_  must  be  taken  into  account  for  intensity  ratio 
techniques . 

The  system  response  was  calibrated  using  a tungsten 
ribbon-filament  standard  of  spectral  radiance.  Figure  11 
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contains  the  results  and  the  system  response  correction 
factors  used  in  the  experimental  work.  Applying  these 
corrections  to  the  raw  data,  values  of  H^lf1  , are  ob- 
tained where  ^ is  a single  valued,  geometry  dependent 
constant . 

There  remains  the  problem  of  relating  X1  to 
^ A.  (X) . For  sources  possessing  cylindrical  symmetry 

/f 

y HK(y)  = 2 f r N(r)r  Jr 

(4.11) 

Figure  12  illustrates  the  geometric  relationships.  Fol- 
lowing the  method  outlined  by  Pearce, ^ a numerical  scheme 
for  solving  equation  (4.11)  for  / (r)  was  developed 

for  the  computer.  The  source  was  assumed  to  be  composed 
of  25  concentric  rings  having  constant  plasma  properties. 
The  code  developed  for  this  work  utilizes  the  actual 
physical  dimension  of  the  source  so  that  absolute  inten- 
sity data  can  also  be  handled  using  the  code.  Because 
of  this  general  utility,  a copy  of  the  program  is  included 
as  an  appendix. 


CHAPTER  V 


EXPERIMENTAL  PROCEDURES 
Experimental  Apparatus 

R.f.  power  was  supplied  by  a Radiomarine  Corporation 
of  America  AN/FRT  15  radio  transmitter.  This  transmitter 
was  tuned  to  operate  in  the  c.w.  mode  at  10  me.  The  power 
was  supplied  to  a tunable  external  circuit  consisting  of  a 
tunable  vacuum  capacitor  and  a 5 turn,  2 inch  diameter 
coil  wound  from  .25  inch  copper  tubing.  This  coil  gen- 
erated the  solenoidal  field  which  excited  the  plasma. 

The  plasma  was  generated  in  a 1.25  inch  diameter,  6 inch 
long,  quartz  tube  which  was  an  appendage  on  a dynamically 
pumped  vacuum  chamber.  The  test  gas  was  leaked  into  the 
vacuum  system  through  a 1/8  inch  Nupro  shutoff  valve 
and  a 1/8  inch  Nupro  very  fine  metering  valve. 

This  system  allowed  the  supply  gas  to  be  shut  off 
completely  for  evacuation  of  the  system  and  at  the  same 
time  provided  fine  adjustment  of  the  operating  pressure 
by  use  of  the  controlled  leak.  The  external  tunable 
circuit  and  the  vacuum  system  were  mounted  together  on 
a K & E instrument  stand.  This  arrangement  is  shown  in 
Figure  13.  An  optical  system  composed  of  a spherical 
front  surfaced  mirror  and  a rotating  flat  front  surfaced 
mirror  mounted  on  an  optical  bench  was  used  to  sweep  a 
full-sized  image  of  the  plasma  across  the  entrance  slit 
of  a modified  Jarrell-Ash  .5  m scanning  monochromator. 
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The  output  of  the  photomultiplier  tube  was  connected  to  a 
Tektronix  547  oscilloscope  using  a type  1A1  amplifier. 

The  rotating  mirror  was  driven  by  a 900  rpm  synchronous 
motor  so  that  the  plasma  (3.175  cm  in  diameter)  was 
scanned  in  0.7372  x 10"3  sec.  This  portion  of  the  ex- 
perimental apparatus  is  shown  in  Figure  14. 

Power  Supply 


The  radio  transmitter  was  operated  at  1/2  maxi- 
mum plate  power.  By  de-tuning  the  external  circuit,  a 
range  of  power  inputs  to  the  plasma  could  be  obtained. 

The  power  to  the  plasma  was  obtained  by  taking  the  dif- 
ference between  the  forward  and  backward  power  as  measured 
by  a Sierra  Philco  bi-directional  power  monitor  and  by 
subtracting  15  watts  which  went  into  heating  the  material 
in  the  neighborhood  of  the  work  coil.  This  power,  lost 
to  the  structure,  was  measured  by  running  the  generator 
while  the  vacuum  system  was  not  operating.  No  plasma 
was  obtained  under  these  conditions  and  the  power  dissi- 
pated in  the  structure  could  be  measured  with  the  power 
meter.  Since  the  load  is  considerably  different  with  no 
plasma  the  dissipated  power  obtained  in  this  manner  is 
only  an  approximation  so  that  the  power  inputs  recorded 
for  the  various  runs  must  be  considered  reasonable  esti- 
mates rather  than  precise  measurements.  Subject  to  these 
limitations,  low  pressure  plasmas  were  generated  with 
input  powers  ranging  from  220  to  740  watts. 

Vacuum  System 


The  vacuum  system  consisted  of  a 2 inch,  Pyrex 
X-section  pumped  by  a Welch  Duo  seal  fore  pump  and  a 2 
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inch  oil  diffusion  pump.  The  diffusion  pump  was  operated 
on  Dow  DC  704  silicon  oil.  Prior  to  a series  of  experi- 
mental runs,  the  system  was  evacuated  to  10-^  mm.  of 
mercury  as  recorded  by  an  ionization  gauge  in  the  vacuum 
line  between  the  diffusion  pump  and  the  chamber.  After 
12  hours  of  operation  at  this  pressure  the  system  was  back 
filled  with  commercial  grade  argon  (99.996%  pure)  through 
the  metering  system  to  the  desired  pressure,  which  was 
measured  using  a thermocouple  gauge  located  in  the  vacuum 
chamber . 

Spectroscopic  Instrumentation 

The  .5  m Jarrell-Ash  scanning  monochromator  had  been 
previously  modified  to  allow  interchange  of  gratings 
through  the  use  of  Kelvin  type  mounts.  Alignment  and 
focusing  of  the  spectrometer  and  the  external  optical 
system  was  carried  out  using  an  Optics  Technology  c.w. 
laser.  The  grating  used  in  this  experiment  was  blazed 
for  7500  S and  had  1180  grooves/mm  which  resulted  in  a 
reciprocal  dispersion  at  the  exit  slit  of  16  S/mm. 

An  EMI  9558B  photomultiplier  was  used  as  the  detector. 

For  all  the  experimental  work,  the  slit  setting  was  100 y. . 
Although  this  did  not  result  in  measurement  of  the 
total  spectral  line  intensity,  it  was  found  that  the  line 
shapes  of  the  lines  involved  were  quite  similar  so  that 
the  ratios  of  the  measured  lines  give  reasonable  values 
for  total  intensity  ratios.  The  alternative  to  the  wide 
slit  would  be  a wavelength  scan  over  the  line  at  each 
spatial  location  and  an  Abel  inversion  for  several  wave- 
length locations  near  the  central  frequency  of  the  line. 
Errors  in  the  magnitude  of  the  intensity  ratios,  due  to 
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neglecting  the  local  line  shapes,  are  expected  to  be  on 
the  order  of  5%.  In  view  of  the  approximate  nature  of 
the  model  used,  this  error  is  not  expected  to  be  a sig- 
nificant  limitation  on  the  accuracy  of  the  temperatures. 

The  signal  output  cable  from  the  photomultiplier  was 
shortened  to  8 inches  and  double  shielded  in  an  attempt 
to  reduce  the  r.f.  noise  which  resulted  from  radiation 
from  the  work  coil.  In  addition,  the  high  voltage  lead 
to  the  photomultiplier  tube  and  the  photomultiplier  tube 
itself  were  shielded  with  copper.  At  high  powers  the  r.f. 
still  imposed  a 10  me  harmonic  on  the  signal  which  ap- 
peared to  be  independent  of  the  gain  on  the  amplifier. 

It  was  assumed  that  this  was  caused  by  the  10  me  signal 
getting  to  the  deflection  plates  in  the  oscilloscope. 

This  noise  was  tolerable  however  and  measurements  could 
be  made.  The  time  constant  of  the  external  circuit  was 
not  quite  small  enough  to  allow  absolute  following  of 
the  signal.  To  minimize  this  effect,  the  output  data 
were  reflected  about  the  centerline  and  an  average  curve 
was  used  in  the  data  analysis. 

The  intensity  distribution  across  the  plasma  torch 
in  a small  wavelength  interval  was  measured  by  sweeping 
a full  scale  image  of  the  torch  across  the  entrance  slit 
of  the  spectrometer  using  the  rotating  mirror.  Alternate 
sweeps  were  used  to  trigger  the  oscilloscope  and  the  sweep 
rate  of  the  delayed  sweep  was  adjusted  to  give  as  large  a 
data  presentation  as  possible. 

Data  Reduction 


The  intensity  distributions  of  four  spectral  emission 
lines  of  argon  were  obtained  for  three  different  input 
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power  levels  at  pressures  of  50  and  100^1 1 . The  four 
lines  used  in  this  study  consisted  of  the  7503.87  X and 
7514.65  X excited  atom  lines  and  the  4806.90  X and 
4347.11  X excited  ion  lines. 

In  addition,  the  intensity  distributions  of  the 
4806.90  X and  7503.87  X argon  lines  were  obtained  at  a 
pressure  of  1000^/  and  an  input  power  of  640  watts. 

Data  obtained  on  Polaroid  photographs  were  response 
corrected  and  plotted  on  graph  paper.  Typical  response 
corrected  data  are  plotted  in  Figures  15,  16,  17  and  18. 

To  get  radial  intensity  distributions,  these  data 
were  used  for  the  numerical  inversions  of  Abel's  equation 
along  with  the  code  found  in  the  appendix.  Four  typical 
radial  intensity  distributions  are  shown  in  Figures 
19,  20,  21  and  22  which  correspond  to  the  previously 
mentioned  response  corrected  data.  The  local  values 
of  intensity  ratios  at  25  radial  locations  for  the  line 

pairS  I7514/  J4347  and  l150Z/  *4806  were  used  to  obtain 
the  temperature  distributions. 


CHAPTER  VI 


RESULTS  AND  CONCLUSIONS 

The  temperature  distributions  in  the  discharge  as 

determined  from  the  line  ratios  I„^„ ./  T and 

7514  4347 

I7503//  *4806  are  Plotted  in  Figures  23  through  29.  All 
of  the  temperature  distributions  have  essentially  the 
same  scatter  as  well  as  the  same  general  shape.  The 
distributions  determined  using  the  ratio  I4g06 

all  exhibit  an  inflection  point  approximately  midway 
between  the  axis  and  the  wall.  This  was  probably  due  to 
error  introduced  in  the  inversion  of  Abel's  equation 
which  causes  some  scatter  near  sharp  peaks  in  the  in- 
verted data. 

The  temperature  distributions  obtained  exhibited 
maximum  temperatures  on  the  axis  of  symmetry  of  the  dis- 
charge. The  magnitude  of  this  on-axis-peak  is  in  ex- 
ceptionally good  agreement  with  the  asymmetric  double  probe 

14 

data  of  Keefer,  et  al.  Averaged,  on-axis  temperatures 
are  presented  in  Table  1.  It  can  be  seen  on  Figure  23 
that  the  temperature  drops  from  this  on-axis  maximum 
to  approximately  85%  of  the  axial  value  at  a distance  of 
2/3  of  the  radius.  No  behavior  of  this  nature  has  been 
observed  with  the  probes.  Referring  to  Figure  30  it  may 
be  seen  that  the  temperature  is  a rather  weak  function 
of  the  input  power  in  the  power  range  of  220  to  740 
watts.  This  particular  behavior  has  not  been  observed 
with  probe  techniques.  It  would  seem  that  the  spectro- 
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metric  technique  is  more  sensitive  to  electron  temperature 
variations  than  the  probe  techniques.  In  view  of  the 
approximate  nature  of  the  plasma  model  used,  the  absolute 
magnitude  of  the  measured  temperatures  is  subject  to  some 
question.  However,  the  general  agreement  between  the 
probe  and  spectrometric  techniques  tends  to  verify  the 
model.  Presumably,  errors  in  calculating  the  excited 
state  populations  in  the  successive  ionization  stages 
tend  to  cancel  when  a ratio  is  taken. 

A calculation  of  the  ratio  of  number  density  of 

Arlll  to  Aril  using  equation  (2.21),  an  assumed  electron 

12  -3 

number  density  of  10  cm  , and  a temperature  of  1.5 
ev.  yields  a number  of  the  order  10^.  Assuming  the 
plasma  to  be  quasi-neutral  and  neglecting  double  ioni- 
zation Dalton's  law  may  be  applied  to  the  plasma.  For 
an  assumed  neutral  gas  and  in  temperature  of  0.2  ev. 

the  electron  number  density  is  approximately  5 x 1013 
-3  . 

cm  for  a total  pressure  of  100  microns.  In  view  of 
the  approximations  made  in  applying  Dalton’s  law  and  the 
approximate  nature  of  the  plasma  model  this  must  be  con- 
sidered an  order  of  magnitude  estimate  for  the  axial 
value  of  the  number  density.  No  attempt  was  made  to 
converge  on  a self-consistent  temperature  and  electron 
number  density  since  the  measured  temperatures  are  only 
slightly  influenced  by  the  electron  number  density 
and  absolute  magnitude  of  the  number  density  would  re- 
main subject  to  the  assumed  ion  and  neutral  temperature. 

The  electron  number  density  could  also  be  inferred 
from  the  absolute  intensity  of  an  ion  emission  line. 
However,  this  requires  an  accurate  value  for  the  popu- 
lation density  of  the  higher  excited  state,  which  the 


32 


approximate  model  used  probably  does  not  provide. 

In  the  type  of  discharge  studied,  it  is  probable 
that  diffusive  loss  of  ions  to  the  wall  is  a significant 
loss  mechanism.  At  the  present  time  there  does  not 
appear  to  be  an  approximate  semi-corona  model  which  will 
suit  such  a loss  mechanism.  If  a solution  of  the  rate 
equations  for  argon  were  possible  this  loss  mechanism 
could  be  accounted  for.  This  seems  to  be  impossible  in 
the  near  future.  It  is  suggested  that  future  work  on  the 
electrodeless  discharge  should  utilize  hydrogen  or  helium. 
The  theory  for  hydrogenic  ions  could  then  be  applied, 
including  the  diffusive  loss  to  the  walls. 
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Figure  1 Temperature  Distribution. 
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Figure  2 Domains  of  Validity  for  an 
Optically  Thin  Argon  Plasma  . 
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Figure  3 Ratio  of  Neutral  Atoms  to 
Singly  Ionized  Ions  for  Argon. 
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Figure  4 Saha  Equation  Correction 
Factor  for  Argon  Atoms. 
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Figure  5 Collision  Limit  for  Argon  Atoms. 
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Figure  6 Collision  Limit  for  Argon  Ions. 


39 


T (ev.) 
e 

Figure  7 Ratio  of  Excited  Neutrals  to  Excited  Ions. 
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Figure  8 Emission  Line  Intensity  Ratio. 
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Figure  9 Law  of  Irradiance. 
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Figure  10  Experimental  Verification  of  the 
Law  of  Irradiance. 
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Figure  11  Spectral  Intensity  Calibration. 
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Figure  12  Geometry  of  the  Source. 
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Figure  15  Integrated  Intensity  Profile  4806  8. 


Voltage 
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Figure  16  Integrated  Intensity  Profile  4347 
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Figure  17  Integrated  Intensity  Profile  7503 
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Figure  18  Integrated  Intensity  Profile  7514 
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Figure  19  Radial  Intensity  Distribution  4806 
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Figure  20  Radial  Intensity  Distribution  4347 
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Figure  21  Radial  Intensity  Distribution  7503 
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Figure  22  Radial  Intensity  Distribution  7514 
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Figure  23  Radial  Temperature  Distribution  Run  13. 
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Figure  24  Radial  Temperature  Distribution  Run  14. 
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Figure  25  Radial  Temperature  Distribution  Run  15. 


GO 


Figure  26  Radial  Temperature  Distribution  Run  16. 
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Figure  27  Radial  Temperature  Distribution  Run  17. 
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Figure  28  Radial  Temperature  Distribution  Run  18. 
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Figure  29  Radial  Temperature  Distribution  Run  19. 
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Figure  30 
Power 


Peak  Temperature  as  a Function  of  Input 


Table  1 


Results 


Run 

No. 

Temperature  (ev) 

Pressure 

Input  Power 
Watts 

*7514 

■*■7503 

Microns  of  Hg 

I4347 

*4806 

13 

1.55 

1.53 

50 

640 

14 

1.36 

1.43 

50 

380 

15 

1.21 

1.26 

50 

220 

16 

1.31 

1.33 

100 

220 

17 

1.47 

1.50 

100 

440 

18 

1.56 

1.53 

100 

740 

19 

1.48 

1000 

640 

65 


APPENDIX 


NUMERICAL  INVERSION  OF  ABEL'S  EQUATION 

The  data  input  for  the  following  Fortran  IV  program 
consists  of  one  card,  using  format  1,  with  the  length  of 
the  source  radius  in  arbitrary  units,  followed  by  as  many 
sets  of  4 cards,  using  format  3,  as  there  are  data  to  be 
processed.  These  4 card  sets  contain  the  integrated 
intensity  data  for  25  radial  locations. 


C PROGRAM  TO  CALCULATE  DISTANCE  FROM  THE  OBJECT  PLANE 

SOURCE  DIVIDED  INTO  25  CONCENTRIC  RINGS 

1 FORMAT  (IF10.5) 

2 FORMAT  (3X, 26HDISTANCE  FROM  OBJECT  PLANE) 

3 FORMAT  (7F10.5) 

4 FORMAT  (3F10.6) 

5 FORMAT  ( 5X,  23HDEFOCUSING  COEFFICIENTS,  3X,  (3F10.6)  ) 

6 FORMAT  (10F10.5) 

12  FORMAT  (IX, 30HDEFOCUSING  COEFFICIENTS  S(K,J)  ) 

DIMENSION  A(26, 26) , RA(25),  R(25) , G(25, 25) , 

E(25,25) , S (25, 25) 

DIMENSION  B(26, 25) , Y(25) , X(26) , Z(25),  V(26,26) 

139  READ  (5,1)  RAD 
DO  7 K = 1,25 
DO  7 J = 1,K 

RA(K)  = (RAD/25.0) *FLOAT(K) 

R(J)  = (RAD/25 . 0) *FL0AT  (J)  - ( (RAD/25 . 0) *2 . 0 - 
(RAD/25.0)  ) *0 . 5 

A(K,J)  = (RA (K) **2.0  - R(J)**2.0)**0. 5 

7 CONTINUE 
WRITE  (6,2) 

WRITE  (6,3)  ( (A (K, J) , J=1 , K) , K=1 ,25) 

C PROGRAM  TO  CALCULATE  LENGTH  COEFFICIENTS  FOR 

ABELS  EQUATION 

15  FORMAT ( IX, 19HLENGTH  COEFFICIENTS) 

DO  120  K = 2,26 
L = K -1 
DO  120  J = 1 ,L 
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120  V(K,J)  = A (L , J) 

DO  125  J = 1,25 
125  V(J,J)  = 0.0 
DO  130  K = 1,25 
DO  130  J = 1 , K 
L = K + 1 

130  B(K, J)  = (V(L , J)  - V(K , J)  )*  2.0 
WRITE  (6,15) 

WRITE  (6,6)  ( (B(K, J) , J = 1,K),K  = 1,25) 

C INVERSION  OF  ABELS  EQUATION  - 25  RINGS,  Y - 

EXPERIMENTAL  DATA 

13  FORMAT  (IX,  15HINPUT  INTENSITY) 

14  FORMAT ( IX,  16HOUTPUT  INTENSITY) 

131  READ  (5,3)  (Y(I),I  = 1,25) 

DO  150  K = 1,25 
I = 26  - K 
X(26)  = 0.0 
J = I 

Z(I)  = 0.0 

140  Z(I)  = Z(I)  + X ( J+l) *B ( J+l , I ) 

IF  (25-J)150,150,145 
145  J = J+l 

GO  TO  140 

150  X( I ) = ( Y( I ) - Z ( I ) )/(B( I , I ) ) 

CONTINUE 
WRITE  (6,13) 

WRITE  (6,6)  (Y ( I ) , 1=1,25) 

WRITE  (6,14) 

WRITE  (6,6)  (X( I) , I = 1,25) 

GO  TO  131 
END 
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